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Abstract 

We calculate the quasiparticle defect states and charge transition levels of oxygen vacancies in 
monoclinic hafnia. The charge transition levels, although they are thermodynamic quantities, can 
be critically dependent on the band gap owing to localized defect states. These quasiparticle defect 
level effects are treated using the first principle GW approximation to the self energy. We show that 
the quality and reliability of the results may be evaluated by calculating the same transition level 
via two physical paths and that it is important to include the necessary electrostatic corrections 
in a supercell calculation. Contrary to many previous reports, the oxygen vacancies in monoclinic 
hafnia are found to be a positive U center, where U is the defect electron addition energy. We 
identify a physical partitioning of U in terms of an electronic and structural relaxation part. 

PACS numbers: 73.40.Qv, 85.30.Tv, 61.72.jd 
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Hafnia has recently received much attention because of its many applications, in particular 
as high-dielectric gate material replacing silica in microelectronic devices. However, devices 
based on hafnia suffer from several problems such as voltage threshold instabilities l| and 
flat band voltage shifts [2]. These problems are believed to be due to a high density of 
defects in the material in particular oxygen vacancies are believed to play an important role 
as electron traps. 

There have been several theoretical studies on the structural and electronic properties 
of oxygen vacancies in monoclinic hafnia. In early studies, formation energies as well as 
defect levels were calculated within density functional theory (DFT) using the local density 
approximation j^] and the generalized gradient approximation (GGA) [4]. The onsite energy, 
U, for adding an additional electron to the defect was calculated and the vacancies were 
foond to be negative U Q centers within these approxrn.at.ons. In these stnd.es, however, 
the defect levels in the gap could not be determined unambiguously, owing to the well-known 
problem of underestimation of band gaps using Kohn-Sham eigenvalues 6j. Later studies 
used hybrid functionals to calculate the defect levels of oxygen vacancies in hafnia. 
These functionals, which were constructed to fix the band gap underestimation problem, 
found that the defect was no longer a negative U center [8(. Recently, there has been a 



higher level theory study 



10( using a combined DFT and GW approach on these defects. 



These authors found a negative U behavior for the oxygen vacancy. The GW part of the 



study in Ref. [10[ was, however restricted to a 24-atom super cell. 



In this letter, we report a new study of the quasiparticle excitations and charge transition 
levels of oxygen vacancies in monoclinic hafnia. To insure the accuracy of the results, we 
computed the charge transition levels via two physical paths and took into account an 
important electrostatic correction that was not included in the previous studies. The defects 
are found to be positive U centers. Further, we develop a methodology to understand the 



defect charging energy. We use a co,„b,ned approach 11] based on DFT and GW for.na.ism 
[6J . This formalism corrects for the error incurred in calculating formation energy and charge 
transition levels within standard DFT. The calculations were done using large super cells 
with 96 atoms. Such large super cells are necessary to minimize any spurious defect-defect 
interactions from overlap of the defect state wave functions in neighboring super cells. 

The formation energy of a defect in charge state q and at arbitrary ionic coordinates R 



and chemical potential /i, E^(R)[/i], can be expressed as 

Ej(R)[^]=E q (R)-E nrf +^q 



(1) 



where E q (R) is the total energy of the system in charge state q and ionic coordinates R, 
and E re f is the energy of a reference system with the same number of atoms as the charged 
system. We note that R is an arbitrary configuration which needs not be the equilibrium 
configuration of the charge state q which we denote as R q . Charge transition level, £ q / q_1 , is 
defined as the value of the chemical potential at which the charge state of the defect changes 
from q to q-1 (q/q-1). Conventionally, one defines the charge transition level from the 
valence band maximum, E v . It is defined as the chemical potential fi at which the formation 
energies of the q and q-1 defects are equal and can be written in terms of formation energies 
as: 

e^- 1 = Eq. 1 (Rq.i)[/x = E v ] - E{(R q )[/x = E v ] (2) 

Within standard DFT, charge transition level is determined by calculating each of the for- 
mation energies in Eq. [2] in their respective equilibrium configurations as accurately as 
possible. But, because of the band gap problem and self-interaction terms within standard 
DFT methods, significant errors may be introduced. However, within the combined DFT 
and GW formalism, a charge transition level is written as: 

e^- 1 HE^fVO - E^(R q )] + [E^(R q ) - E^(R q )] 

=E re lax + Eqp (3) 

by adding and subtracting the term Eq_ 1 (R q ) |rjj. (All the formation energies in the above 
expression as well as throughout the rest of the letter are evaluated with fi = E v .) This 
reformulation allows us to combine terms to eliminate most of the errors mentioned above. 
Figure [1] shows schematically the procedure for calculation of charge transition levels within 
the DFT+GW formalism. The first bracketed term on the right hand side of the Eq. [3] is a 
relaxation energy E re iax (red line in Figure H]) and the second bracketed term is a quasiparticle 
excitation energy Eqp (blue line in Figure H]). For the excitation energy, one uses the GW 



formalism jfj]; while for the relaxation energy, one uses DFT. This ensures an accurate 
calculation of the appropriate physical quantities and hence the charge transition level and 
relative formation energies. 
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FIG. 1. Formation energies evaluated with fx = E v vs. generalized coordinate, illustrating the 
terms in the DFT+GW formalism for the charge transition level e q//q_1 

It is worth noting that the charge transition level is a thermodynamic quantity and does 
not depend on the path in the formation energy-generalized coordinate space one takes to 
calculate it. In other words, the value of charge transition level remains unaffected when one 
adds and subtracts any formation energy to it. In particular, we can alternatively choose 
to add and subtract E^(R q _i) in Eq. 121 This would correspond to another path (the green 
line) in Figure [TJ Calculating the charge transition level via multiple paths not only serves 
as a check for our calculations, but also gives some idea about the accuracy of the method. 

For the DFT part of our calculation we used an ab initio pseudopotential plane-wave 



method, as implemented in Quantum ESPRESSO |l2j, with PBE 13| exchange correla- 



tion functional. We used non-local pseudopotentials constructed using the Troullier Martins 
scheme with valence configurations 5s 2 5p 6 5d 2 6s 2 and 2s 2 2p 4 for Hf and O, respectively. 
The electronic wave functions were expanded in plane waves with cutoff energy of 250 Ry. 
The k-point sampling was restricted to the T point in view of the large super cell used. 
Our calculated lattice parameters for monoclinic hafnia are in excellent agreement with 



experiment [15| as well as with previous calculations |4J. Using the PBE eigenvalues and 
eigenfunctions, the quasiparticle energies were calculated within the GnWp approximation 
[6] to the electron self- energy as implemented in the BerkeleyGW package (l6j. The static 
dielectric matrix was calculated with a 10 Ry energy cutoff in a plane wave basis and ex- 
tended to finite frequencies within the generalized plasmon pole model pj. The band gap for 
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rx V~5 

FIG. 2. (a) Schematic diagram showing both three and four fold coordinated oxygen vacancy 
induced defect levels in the band gap of hafnia for various charge states, (b) Periodic images of 
charged defects leading to a spurious electrostatic potential on the defect, (c) Effect of this spurious 
potential on eigenvalues. The plot shows cell size vs PBE eigenvalue with respect to the top of the 
valence band for vacancy induced defect states in Vo 3 - The dashed line is the linear extrapolation 
to infinite super cell size. 



bulk monoclinic ha: 
studies 5.45-5.9 



nia is calculated to be 6.0 eV, which is in_good agreement with previous 



Il7j ] as well as experiments 5.7-5.9 [l8|, Il9 |. 

There are two distinct types of oxygen vacancies (with different coordination) in mon- 
oclinic hafnia - 3-fold coordinated (Vo 3 ) and 4-fold coordinated (Vo 4 ). We performed 
calculations for charge states q = 0, 1, 2 for both kind of vacancies, i.e. zero, one and two 
missing electrons from the vacancy. 

Figure [2] (a) shows a schematic of the oxygen vacancy induced quasiparticle defect levels 
in the band gap of hafnia. The defect levels for all the charge states lie deep in the gap. Table 
fl] shows the relaxation energies (E re i ax ), quasiparticle energies (Eqp) and charge transition 
levels (e q / q_1 ) calculated within DFT+GW approach. Despite the large differences in the 
quasiparticle energies required for the calculation in the two paths, the final charge transition 
levels via the two paths are always within 0.2 eV of one another. This gives us an error 
estimate of ±0.1 eV for the DFT+GW method for calculation of charge transition levels. 

It is worth noting that quasiparticle energies in Table [I] include an electrostatic correction 
owing to the super cell geometry used. Within the standard DFT-only methodology of 
calculating charge transition levels, one corrects the total energies of charged defects for 
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-0.64 0.76 


-0.57 0.67 


-0.75 0.80 


-0.55 0.65 


E QP 


3.30 1.69 


3.93 2.88 


3.04 1.33 


3.46 2.50 


e q/q-l 


2.66 2.45 


3.36 3.55 


2.29 2.13 


2.91 3.15 


Avg. e^" 1 


2.56 


3.46 


2.21 


3.03 



TABLE I. Table showing the two components contributing to the charge transition levels for Vo 3 
and Vo 4 calculated within the DFT+GW methodology. PI (P2) refers to Path 1 (2). Row 1 shows 
the relaxation energy (E re i ax ), Row 2 the quasiparticle energy (Eqp), Row 3 the charge transition 
level (e q / q_1 ) and Row 4 the value for e^^~ l averaged over both paths. All values are in eV. 



unphysical electrostatic terms from the charge on the image defects using Makov-Payne 
201 ] like corrections. Within the DFT+GW method, the electrostatic error in the DFT 

eigenvalues in general needs to be accounted for. The origin of these electrostatic errors in 
;he DFT eigenvalues is the spurious electrostatic potential from charged neighboring defects 



2l| as shown schematically in Figure |2] (b). To quantify these errors, we plotted the position 



of the Kohn-Sham defect level with respect to the valence band maximum for various super 
cell sizes as shown in Figured (c). These calculations were done using the SIESTA code 
22 ] with increasing super cell sizes - the largest of which contained 1499 atoms + defect. To 
get the infinite cell size value, we linearly extrapolated from the largest two super cell sizes. 
Figure[2](c) shows that even at 1500 atoms, the errors are large with respect to the converged 
answer. Table [Til shows the effect of electrostatic corrections on quasiparticle defect levels 
along the two paths for Vo 3 - The two paths have different electrostatic corrections to the 
quasiparticle energies and that these corrections are substantial. Accounting for these errors 
is important to get a reliable value of charge transition level. 

Our quasiparticle defect levels and charge transition levels are in disagreement with pre- 
vious DFT+GW calculations jlQj ]. This disagreement is likely a consequence of their choice 
of a small super cell and, more importantly, neglect of electrostatic corrections. Reference 
Q's value of e +2 / +1 for Vo 3 (V J is 4.00 eV (3.22 eV) and of e +l l° is 3.10 eV (2.43 eV); 
these values are close to our corresponding uncorrected values. Our calculations also dis- 



+2/+1 


+1/0 


PI P2 


PI P2 


E QP (Uncorrected) 5.14 2.66 


4.90 2.88 


E QP (Corrected) 3.30 1.69 


3.93 2.88 


AEqp 1.84 0.97 


0.97 0.00 



TABLE II. Table showing effect of spurious charged super cell electrostatic corrections on the 
quasiparticle defect levels for Vo 3 (referred to the top of the valence band). PI (P2) refers to Path 
1 (2). Row 1 shows the uncorrected quasiparticle energy, Row 2 the corrected quasiparticle energy, 
Row 3 the correction. All values are in eV. 
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FIG. 3. Relative formation energy vs. chemical potential for charged oxygen vacancies in mono- 
clinic hafnia. The shaded region denoted the Si band gap. Thee formation energy of Vq 4 is lower 
than V^ 3 by 0.14 eV. 



of 3.7 eV for s +2 l +1 in Vq 3 and 4.0 eV for 



agree with the hybrid functional value 
between Vq 3 and Vq 4 . (Our corresponding values are 2.56 eV and 3.32 eV respectively.) It 
is also worth noting that the Kohn-Sham eigenvalues calculated with hybrid functionals also 
disagree with our calculations. Also, in hybrid calculations, even though total energies were 
corrected for electrostatic errors, the eigenvalues were not. While this is not necessary for 
the charge transition level calculated as a total energy difference, it can lead to erroneous 
results whenever eigenvalues are directly compared to single-particle excitation experiment. 



Figure |3] shows the relative formation energy of various charge states for the oxygen 
vacancies as a function of the chemical potential in monoclinic hafnia based on the values in 
Table [H The formation energy in Figure [3] is plotted with respect to the formation energy of 
V& 3 . As evident from Figure El V^, and V^ 4 are the most stable defects in the system 
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depending on the value of /i. The experimentally relevant charge transition levels, e +2 l +1 
and would be 2.56 eV and 3.32 eV respectively. Not shown in the figure is the stability 
for VoVVq 2 defects. For larger values of the chemical potential, the VqVVq 2 defects may 
become more stable. Also, shown in Figure [3] (in shaded grey), is the band gap and the 
expected band offset for Si (~3 eV) wffen p.aeed next to hafnia Q. For p-doped Si next 
to hafnia, the system is expected to have Vq* vacancies. 

It is also noted from Figure |3] that oxygen vacancies are positive U centers. U is defined 
as the energy of the reaction: 2V +1 — > V 2+ + V°. In terms of the charge transition levels, U 
for V +1 can be written as: 



Our calculated values of U for both vacancies (Vq* and Vq*) are given in Table HTT1 Further, 
U can be broken into two parts - an electronic part (U e i ec ) and a relaxation part (U rc i ax ) as 
defined in the curly brackets in the Eq. 5 below: 



This partitioning of U is instructive, because physically, U e i ec represents the quasiparticle 
gap of the system in the +1 charge state keeping the structure fixed and U re i ax represents 
sum of structural relaxation energies. Physically, U e i ec > and U re i ax < 0. Table II III shows 
our calculated values of U re i ax and U e i ec - The reason for large relaxation energy is that, in 
the +1 and +2 charge state, the atoms nearest to the vacancy relax by up to 5-10% of their 
)ond lengths. This is in agreement with previous studies of relaxation around the vacancies 



In conclusion, we have reported quasiparticle energies and charge transition levels of 
oxygen vacancies in monoclinic hafnia. We find that Vq 2 , Vo* and Vq 4 are the most stable 
oxygen vacancies in the system as the Fermi level spans the band gap. By calculating the 
charge transition levels via two paths in configuration-space we gain insight to the charge 
defect stability and highlighted the importance of electrostatic corrections in super cell defect 
calculations. Further, we developed an intuitive partitioning of the defect charging energy U 



U = E^ 2 (R +2 ) + e£(R ) - 2E f +1 (R +1 ) 



(4) 



U ={E f +2 (R +1 ) + E f (R +1 ) - 2Ei 1 (R +1 )}+ 

{[E f +2 (R +2 ) - E f +2 (R +l )} + [EftRo) - EjCR+0]} 



(5) 




Urn- 
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U e lec U re lax Ut t a ] 

2.24 -1.33 0.90 
V+l 2.13 -1.35 0.81 

TABLE III. U for V+J and calculated using the charge transition levels from Table [H Also 
shown are the contributions to U from electronic and relaxation components. 

into a quasiparticle gap and sum of relaxation energies. Contrary to some previous studies, 
which found negative U or U close to zero, the oxygen vacancies were found to be positive 
U centers. 
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